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Background: Recent studies have suggested that several ovarian cancer risk factors differ by parity status, but these findings have 
not been confirmed. We evaluated whether known risk factors of ovarian cancer differ between nulliparous and parous women 
using data from two large prospective cohorts. 

Methods: Data from the National Institutes of Health-AARP Diet and Health Study and the Prostate, Lung, Colorectal and Ovarian 
Cancer Screening Trial were combined for this analysis. Cox regression models were used to estimate associations with ovarian 
cancer risk. Risk heterogeneity by parity status was assessed using likelihood-ratio tests. 

Results: Among the 125437 women included in the analysis, there were 16589 (13%) nulliparous women and 108848 (87%) parous 
women. Of the 623 women diagnosed with invasive epithelial ovarian cancer, 102 (16%) were nulliparous and 521 (84%) were 
parous. While parity reduced ovarian cancer risk, no differences were found for other risk factors by parity. Among ever users of 
hormone therapy, body mass index suggestively increased the risk of ovarian cancer by 1.5-fold in nulliparous but not parous 
women (P-heterogeneity = 0.08). 

Conclusion: While nulliparous women have higher ovarian cancer risk than parous women, our findings suggest that the relative 
effects of most other risk factors do not differ by parity. 



Parity is a well-established protective factor for ovarian cancer 
(Hankinson and Danforth, 2006). In a collaborative analysis of 
12 case-control studies, women with at least one full-term 
pregnancy were at a significantly lower risk compared with 
nulliparous women (odds ratio (OR) = 0.60, 95% confidence 
interval (CI): 0.47-0.76), and each additional full-term pregnancy 
lowered the risk further (P trend < 0.001) (Whittemore et al, 1992). 
Similar results have been found in cohort studies (Adami et al, 
1994). 

Recently, it has been suggested that associations between several 
ovarian cancer risk factors, such as body mass index (BMI) or the 
use of oral contraceptives (OCs), may differ by parity status (Ness 
et al, 2001; Greer et al, 2006). Greer et al (2006) found that a high 
BMI significantly increased the risk of ovarian cancer in 
nulliparous women (OR = 2.53, 95% CI: 1.39-4.61 comparing 



the top to the bottom quartile of the BMI distribution), but not in 
parous women (OR=0.96, 95% CI: 0.70-1.31 comparing the top 
to the bottom quartile of the BMI distribution). In the same 
population, use of oral contraceptive (OC) did not reduce the risk 
of ovarian cancer in nulliparous women (OR =0.9, 95% CI: 
0.5-1.7) but was suggestively associated with a reduced risk in 
parous women (OR = 0.6, 95% CI: 0.3-1.1 in women with one 
birth; OR = 0.6; 95% CI: 0.4-1.0 in women with two births; 
OR=0.7, 95% CI: 0.5-1.0 in women with three or more births) 
(Ness et al, 2001). The authors did not explore whether other 
ovarian cancer risk factors differ by parity. These results were 
based on a retrospective sample, which could have affected their 
BMI results. In addition, the study population included pre- and 
post-menopausal women, with 40% of the cases being younger 
than 50 years of age (Ness et al, 2000), and associations between 
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BMI or OC use and ovarian cancer risk may differ between pre- 
and post-menopausal women (Tung et al, 2005; Moorman et al, 
2008; Olsen et al, 2013). A further limitation of the studies by Ness 
et al (2000, 2001) is the inclusion of borderline and invasive 
ovarian cancer cases. It has been suggested that BMI and OC use 
could be more strongly associated with borderline tumours 
(Modugno et al, 2001; Collaborative Group on Epidemiological 
Studies of Ovarian Cancer, 2012; Olsen et al, 2013). Therefore, 
these results need confirmation in large, prospective data sets, with 
information on menopausal status and use of hormone therapy 
(HT) among postmenopausal women. A different study reported 
an interaction effect between family history and parity on the risk 
of ovarian cancer (Vachon et al, 2002). 

To our knowledge, a systematic investigation of whether parity 
modifies the association between other ovarian cancer risk factors 
and ovarian cancer risk has not been done. Understanding the 
interplay between parity and other ovarian cancer risk factors 
could help to better understand the biological mechanism in the 
development of ovarian cancer. In the current analyses, we, 
therefore, evaluated whether the associations of several known risk 
factors for ovarian cancer differ between nulliparous and parous 
postmenopausal women using combined data from two large 
prospective cohorts. 



MATERIALS AND METHODS 



Study population. Data from two prospective cohorts, the 
National Institutes of Health-AARP Diet and Health Study 
(NIH-AARP) (Schatzkin et al, 2001) and the Prostate, Lung, 
Colorectal, and Ovarian Cancer Screening Trial (PLCO) (Prorok 
et al, 2000) were combined for this analysis (Figure 1). Both studies 
have been previously described (Prorok et al, 2000; Schatzkin et al, 



2001). Briefly, the NIH-AARP study was established in 1995-1996 
when over 500 000 AARP members, aged 50-71 in six states 
(California, Florida, Louisiana, New Jersey, North Carolina and 
Pennsylvania) and two metropolitan areas (Atlanta, GA, and 
Detroit, MI) satisfactorily completed a baseline questionnaire 
(Schatzkin et al, 2001). A second questionnaire was sent in 1996- 
1997 to obtain additional risk factor information, including family 
history of breast and/or ovarian cancer. Information from both 
questionnaires was used in this analysis (baseline questionnaire 
collected in 1995-1996 and second questionnaire collected in 
1996-1997). The last follow-up date for this study was 31 
December 2006. The NIH-AARP study was approved by the 
Special Studies Institutional Review Board of the National Cancer 
Institute. The PLCO study is a multi-centre, randomised trial 
designed to evaluate screening modalities for prostate, lung, colon, 
and ovarian cancers (Prorok et al, 2000). Participants were 
recruited between 1992 and 2001 and they were followed for up 
to 13 years. For this analysis, the last follow-up time for the PLCO 
was 31 December, 2006. At enrolment, women aged 55-74 years 
old, without a history of ovarian cancer completed a baseline 
questionnaire. Exposure data from this questionnaire, collected 
from 1992 to 2001, was used in the analysis. Only information 
from women in the screening arm of the PLCO study was available 
when these analyses were conducted. Informed consent for the 
PLCO was obtained from all women and each screening centre 
obtained approval by the local institutional review boards. In both 
studies, participants completed self- administered questionnaires or 
telephone-based interviews with information regarding demo- 
graphic and lifestyle characteristics, as well as reproductive and 
medical history. Data from the both studies were harmonised into 
a single data set (Supplementary Note). Analyses were restricted to 
women who reported at age at menopause at the beginning of 
either study. 



138 056 Women in NIH-AARP study 
(Second baseline questionnaire) 



1649 Women that used a proxy 
respondent were excluded 




35 552 Women in PLCO screening arm 



7824 Women with missing age at menopause 
were excluded 



1 64 1 35 Women in combined data set 



37 541 Women with both ovaries removed or 
unknown status were excluded 



126 594 Women in combined data set 



488 Women with a personal history of ovarian 
cancer or prevalent ovarian cancer at study entry 
were excluded 



126 106 Women in combined data set 



669 Women with missing information on parity 
were excluded 




Figure 1. Flow diagram of criteria and numbers of study participants included in the analysis. Of the 125 437 women included in the 
analysis, 30 197 (24.1%) came from the PLCO study and 95240 (75.1%) came from the NIH-AARP study. Reflecting these proportions, of the 
623 ovarian cases, 153 (24.6%) came from the PLCO study and 470 (75.4%) from the NIH-AARP study. 
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Outcome ascertainment. National Institutes of Health-AARP 
participants were followed through the US Postal Service national 
database for address changes and cancers were ascertained by 
linkage to state cancer registries in the original recruitment areas 
plus three catchment areas that were common states of relocation 
(Arizona, Nevada, and Texas) (Schatzkin et al, 2001). Vital status 
was updated from the US Social Security Administration Death 
Master File and the National Death Index Plus (Schatzkin et al, 
2001). The completeness of case ascertainment has been previously 
reported to have a minimum sensitivity of ~ 90% and specificity of 
99.5% compared with identification of cases by cancer registries 
(Michaud et al, 2005). Follow-up in the PLCO was conducted from 
annual study update forms completed by participants or a proxy 
and subsequently confirmed from medical records, death certifi- 
cates, and cancer registries (Prorok et al, 2000). Subjects who did 
not return the questionnaire were contacted by repeat mailing or 
telephone to improve completeness of end point ascertainment. 

The outcome of interest for this analysis was primary invasive 
epithelial ovarian cancer. Using the International Classification of 
Diseases for Oncology, 3rd edition, (ICD-O-3) coding system 
(Fritz et al, 2000), these were cases with C56.9 site codes and 
epithelial morphology (serous: 8441, 8460, 8461; endometrioid: 
8380, 8381, 8560, 8570; mucinous: 8470, 8471, 8480, 8481; clear 
cell: 8310, 8313; and other epithelial: 8010, 8020, 8021, 8050, 8070, 
8120, 8140, 8240, 8246, 8255, 8260,8323, 8440, 8450, 8490, 8562). 

Statistical methods. Cox regression models were used to compute 
hazard ratios (HRs) and 95% CIs as measures of association 
between risk factors and incidence of ovarian cancer. Follow-up 
started at the age at the time of the second questionnaire for 
women in the NIH-AARP study or the age at the baseline 
questionnaire for women in the PLCO and ended at the age of 
diagnosis of invasive epithelial ovarian cancer or censoring at the 
age of death, loss of follow-up, or 31 December 2006, whichever 
occurred first. The baseline hazard function was stratified by study 
(PLCO vs AARP) in all models. 

Women who reported at least one live birth or provided an age 
at first birth were classified as parous and those who reported no 
age at first birth and no live births were considered nulliparous, 
similar to Schonfeld et al (2011). Subjects were excluded from the 
analysis if they did not report age at menopause at the beginning of 
the study, had a bilateral oophorectomy or unknown ovarian 
surgery status, had a personal history of ovarian cancer, prevalent 
ovarian cancer at study enrolment, had missing information on 
parity, or if questionnaire information was obtained via proxy 
respondents (see Figure 1). 

Analyses were done separately for nulliparous and parous 
women from both cohorts combined. Analyses for nulliparous 
women were adjusted for BMI at study entry (continuous), 
duration of use of OC (never or < 1/1—9/ ^ 10 years), duration of 
use of HT (never/ < 10/ > 10 years), first degree family history of 
breast and/or ovarian cancer (no/yes/missing). Analyses for parous 
women were adjusted for the same variables as nulliparous 
women and in addition for the number of live births in categories 
(1/2/3-4/^5). When models were additionally adjusted for 
education, marriage, age at menopause, age at menarche or 
hysterectomy status results did not significantly changed. To assess 
differences in effects by parity status, we compared a Cox model 
that combined the data of the nulliparous and parous women and 
was adjusted for all the above factors and parity (no/yes) to a 
similar model with an additional interaction term between parity 
and the risk factor of interest modelled as a continuous variable. 
Significance of the interaction term was determined by a 
likelihood-ratio test. 

The number of lifetime ovulatory cycles (LOCs) was computed 
using the model by Cramer et al (1995) that estimates LOC based 
on the age at menopause, age at menarche, time being pregnant, 



duration of OC, and the average cycle length. We assumed that the 
average cycle length of all participants was 28 days. Estimated 
LOCs from this model result in similar estimates as to those from 
other published models (George, 2011). Quartiles of LOCs 
were defined based on their distribution in the entire study 
population. 

The proportional hazard assumption was tested based on the 
slope of the Schoenfeld's residuals (Grambsch and Therneau, 
1994). All statistical tests were two-sided and P-values <0.05 were 
considered statistical significant. Analyses were performed by using 
STATA, version 10.1, software (StataCorp LP, College Station, TX, 
USA). 



RESULTS 



A total of 125 437 women were included in the analysis, with 16 589 
(13%) nulliparous women and 108 848 (87%) women who reported 
having had at least one live birth. During the study period, 623 
women were diagnosed with invasive epithelial ovarian cancer; 102 
(16%) were nulliparous and 521 (84%) reported having had at least 
one live birth. The median follow-up time for women in this 
analysis was 10.55 years (interquartile range: 10.52, 10.80 years). 

Characteristics of study participants according to parity status 
are shown in Table 1. Nulliparous women were younger, had 
higher education, were more likely to never have been married and 
were less likely to have used OCs, HT or to have undergone 
hysterectomy or tubal ligation compared with parous women. 
Parous women were at a reduced risk of ovarian cancer 
compared with nulliparous women (HR = 0.79, 95% CI: 0.63, 
0.98) (Table 2), with stronger reduction of the risk with more 
births (P trend < 0.001). 

There were no clear differences in the risk estimates for 
hormonal factors and other characteristics between parous and 
nulliparous women (Table 3). The use of OC reduced the risk of 
ovarian cancer similarly in both nulliparous and parous women 
(P-heterogeneity = 0.98). Similar results were observed when 
analyses were restricted to nulliparous women who reported ever 
having been married (married, widowed, divorced, or separated) or 
living as married. Long duration of use of HT increased the risk in 
both groups of women (Table 3). Compared with women who 
had an estimated number of 470 or more LOCs, both nulliparous 
and parous women with an estimated number of < 330 LOCs had 
lower risk of ovarian cancer. First degree family history of breast 
and/or ovarian cancer increased the risk of ovarian cancer in both 
groups of women (HR= 1.37, 95% CI: 0.81, 2.34 in nulliparous 
women; HR= 1.21, 95% CI: 0.97, 1.51 in parous women) 
(P-heterogeneity = 0.60). BMI was associated with a non-statisti- 
cally significant increased risk of ovarian cancer in nulliparous 
women (HR= 1.36, 95% CI: 0.79, 2.36, comparing 25-29.9 kgm~ 2 
to <25kgm~ 2 ; HR=1.48, 95% CI: 0.61, 3.58 comparing 
30kgm~ 2 or more to <25kg/m 2 ), but was unrelated to ovarian 
cancer risk in parous women (HR = 0.84, 95% CI: 0.66, 1.07, 
comparing 25-29.9 kgm~ 2 with <25kg/m~ 2 ; HR=1.03, 95% 
CI: 0.69, 1.53 comparing 30kgm~ 2 or more with <25kgm~ 2 ) 
(P-heterogeneity = 0.15). BMI was still associated with a non- 
significant increased ovarian cancer risk among ever married 
nulliparous women although the association was attenuated 
(HR= 1.33, 95% CI: 0.51, 3.47 comparing 30kgm~ 2 or more to 
<25kgm~ 2 ) (P-heterogeneity = 0.19). Compared with never 
smokers, nulliparous women who were current smokers were at 
an increased risk of ovarian cancer (HR= 1.40, 95% CI: 0.78, 2.49), 
but the increased was not seen for parous women who were current 
smokers (HR=0.76, 95% CI: 0.55, 1.05) (P-heterogeneity = 0.10). 
Similar results were obtained when analyses were further adjusted 
by education level. 
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Table 1. Distribution of characteristics of subjects at baseline by parity 





Nulliparous 






women 


Parous women 




(N = 16 589) 


(N=108 848) 


Characteristics 


n (%) 


n (%) 


Age (years) 3 


<55 


2273 (13.7) 


8462 (7.8) 


55-<60 


4510 (27.2) 


27 277 (25.1) 


60-<65 


4294 (25.9) 


33 008 (30.3) 


>65 


5512 (33.2) 


40101 (36.8) 


Race/Ethnicity 


Non-hispanic white 


15176 (91.5) 


99 205 (90.7) 


Non-hispanic black 


653 (3.9) 


5336 (4.9) 


Hispanic 


229 (1 .4) 


1719 (1.6) 


Other 


395 (2.4) 


2138 (2.0) 


Missing 


136 (0.8) 


948 (0.9) 


Education 


Less than high school 


2823 (17.0) 


27 280 (25.1) 


High school 


1609 (9.7) 


16 750 (15.4) 


Some college/vocational 


3570 (21.5) 


30893 (28.4) 


College graduate or more 


8285 (49.9) 


31 718 (29.1) 


Missing 


302 (1 .8) 


2207 (2.0) 


BMI (kgrrT 2 ) 3 


<25 


7842 (47.3) 


47 383 (43.5) 


25.0-29.9 


4804 (29.0) 


35 077 (32.2) 


> 30.0 


3517 (21.2) 


23 963 (22.0) 


Missing 


426 (2.6) 


2425 (2.3) 


Marital status 


Married or living as married 


4665 (28.1) 


57 844 (53.1) 


Widowed 


2092 (12.6) 


24249 (22.3) 


Divorced 


2597 (15.7) 


24059 (22.1) 


Separated 


109 (0.7) 


1 563 (1 .4) 


Never married 


7037 (42.4) 


549 (0.5) 


Missing 


89 (0.5) 


584 (0.5) 


Duration of use of oral contraceptives 


Never users or < 1 year 


11 932 (71.9) 


64061 (58.9) 


1-9 years 


3290 (19.8) 


32 730 (30.1) 


S 10 years 


1318 (8.0) 


11 520 (10.6) 


Missing 


49 (0.3) 


537 (0.5) 


Duration of hormone therapy use 


Never users 


9003 (54.3) 


52 336 (48.1) 


< 10 years 


5466 (33.0) 


38 685 (35.5) 


S 10 years 


2064 (12.4) 


17413 (16.0) 


Missing 


56 (0.3) 


414 (0.4) 


Age at menarche 


Early 


7505 (45.2) 


44 409 (40.8) 


Middle 


7278 (43.9) 


49 433 (45.4) 


Late 


1765 (10.6) 


14818 (13.6) 


Missing 


41 (0.3) 


188 (0.2) 


Age at last menstrual cycle (years) 


s:44 


4027 (24.2) 


30 996 (28.3) 


45^9 


4599 (27.6) 


25 583 (23.4) 


50-54 


6772 (40.7) 


41 446 (37.9) 


S55 


1191 (7.2) 


10 823 (9.9) 


Hysterectomy 3 


No 


14123 (85.1) 


77 921 (71.6) 


Yes 


2 410 (14.5) 


30 569 (28.1) 


Missing 


56 (0,3) 


358 (0.3) 



Table 1. (Continued) 


Characteristics 


Nulliparous 

women 
(N=16 589) 
n (%) 


Parous women 
(N = 108 848) 
n (%) 


Tubal ligation 13 


No 
Yes 

Missing 


2524 (92.8) 
1 74 (6.4) 
22 (0.8) 


21 084 (76.7) 
6305 (23.0) 
88 (0.3) 


Endometriosis'" 


No 
Yes 

Missing 


2382 (87.6) 
227 (8.3) 
111 (4.1) 


24 895 (90.6) 
1502 (5.5) 
1080 (3.9) 


Miscarriages/abortions b 


No 
One 

Two or more 
Missing 


2250 (82.7) 
300 (1 1 .0) 
168 (6.2) 
2(0.1) 


17 593 (64.0) 
6383 (23.2) 
3463 (12.6) 
38 (0.1) 


Infertility 1 " 






No 
Yes 

Missing 


1917 (70.5) 
785 (28.9) 
18 (0.7) 


24070 (87.6) 
3324 (12.1) 
83 (0.3) 


First degree family history of breast and/or ovarian cancer 


No 
Yes 

Missing 


10429 (62.9) 
2891 (17.4) 
3269 (19.7) 


77 367 (71.1) 
19729 (18.1) 
11 752 (10.8) 


Smoking status 3 


Never 
Former 
Current 
Missing 


7676 (46.3) 
6344 (38.2) 
2238 (13.5) 
331 (2.0) 


51 083 (46.9) 
41 830 (38.4) 
13 847 (12.7) 
2088 (1 .9) 


Abbreviations: BMI = body mass index; PLCO = Prostate, Lung, Colorectal, and Ovarian 
Cancer Screening Trial. Data might not add up to 100% because of rounding. Age at 
menarche groups for the AARP data set: Early, ^12 years; Middle, 13-14 years; Late, S= 1 5 
years. Age at menarche groups for the PLCO data set: Early, ^ 1 1 years; Middle 1 2-1 3 years; 
Late, >14 years. 
a At enrolment. 

^Only available in the PLCO data set. 

infertility was defined as trying to get pregnant for 12 months or more without success. 



We further explored whether the relationship between BMI and 
parity differed by the number of births among parous women 
(Supplementary Table 1). Among the 108 370 parous women 
with information on number of births, there were 11 706 (10.8%) 
with one birth, 31 451 (29.0%) with two births, 48 397 (44.7%) with 
three or four births, and 16 816 (15.5%) with five or more births. 
Five-hundred seventeen invasive ovarian cancer cases among 
parous women had information on the number of births. Of these, 
75 (14.5%) had one birth, 158 (30.6%) had two births, 225 (43.5%) 
had three to four births and 59 (11.4) had five or more births. 
Comparing nulliparous women with women with one, two, three 
to four or five or more births, there was no clear pattern of 
effect modification (P-heterogeneity = 0.05, 0.63, 0.14, and 0.35, 
respectively). 

As HT use has been suggested to modify the association 
between BMI and ovarian cancer, we additionally stratified the 
BMI analyses by HT use (Table 4). Among never HT users, the 
association between BMI and ovarian cancer was stronger in 
nulliparous women compared with parous women, although the 
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Table 2. Risk of invasive epithelial ovarian cancer in relation to parity 



Risk factor 


Cases (N=623) n (%) 


HR a 


(95% Cl) a 


Parity 


Nulliparous 
Parous 


95 (16.1) 
496 (83.9) 


1.00 

0.79 


ref. 
(0.63, 0.98) 


Number of births 


Nulliparous 

1 

2 

3-4 

>5 


95 (16.1) 
71 (12.1) 
155 (26.3) 
216 (36.7) 
52 (8.8) 


1.00 
1.09 
0.87 
0.75 
0.51 


ref. 
(0.80, 1 .48) 
(0.67, 1.13) 
(0.59, 0.96) 
(0.36, 0.72) 








Ptrend<0.001 


Abbreviations: CI —confidence interval; HR — hazard ratio; PLCO = Prostate, Lung, 
Colorectal, and Ovarian Cancer Screening Trial. Numbers only for women with no-missing 
values for the exposure of interest and corresponding adjusting variables. 
a Adjusted for BMI (continuous), duration of use of oral contraceptives (never or <1/1-9/ 
> 10 years), duration of use of hormone therapy (Never/ <10/> 10 years), first degree family 
history of breast and/or ovarian cancer (no/yes/missing) and stratifying the baseline function 
by study (PLCO vs AARP). 



difference was not statistically significant (P-heterogeneity = 0.76). 
Among ever HT users, nulliparous women with BMI above 
30 kg m ~ 2 had a non-statistically significant higher risk of ovarian 
cancer compared with women with BMI <25kgm~ 2 (HR= 1.50, 
95% CI: 0.57-3.92). This increased risk was not seen for parous 
women (HR = 0.70, 95% CI: 0.44-1.11) (P-heterogeneity = 0.08). 

Finally, we did not observe a difference between histological 
subtypes of epithelial ovarian cancer (serous, mucinous, endome- 
troid, clear cell, and other) and parity status (nulliparous vs parous) 
(/ 2 -test P-value: 0.13) (Supplementary Table 2). 



DISCUSSION 



We assessed whether the relationship between hormone-related 
factors and ovarian cancer risk differed between nulliparous and 
parous women using data from two large prospective cohorts. 
Overall, we did not observe significant differences in risk factor 
associations between nulliparous and parous women. We only 
found a suggestive difference of the associations between BMI and 
ovarian cancer risk in ever users of HT, whereby the estimated 
ovarian cancer risk associated with BMI was greater among 
nulliparous women than parous women. However, because of the 
number of tests performed, this finding could be due to chance. 
These data suggest that the impact of most factors on ovarian 
cancer risk is the same for nulliparous women as for parous 
women, even though ovarian cancer incidence is higher among 
nulliparous women. 

We found a suggestive difference between never vs current 
smokers and risk of ovarian cancer by parity. However, in a recent 
meta-analysis of 51 epidemiological studies, including the NIH- 
AARP study and the PLCO trial, there was no difference in ovarian 
cancer risk between never and current smokers when stratified by 
parity (Collaborative Group on Epidemiological Studies of Ovarian 
Cancer et al, 2012). Given the small numbers of ovarian cancer 
cases who were current smokers in our study, this suggestive 
difference could just be chance. 

Previously, it had been found that BMI significantly increased 
the risk of ovarian cancer by a 2.5-fold in nulliparous women, but 
not in parous women (Greer et al, 2006). We also observed that 
BMI increased the risk of ovarian cancer by nearly 1.5-fold in 
nulliparous women, albeit non-significant, but not in parous 



women. However, in both studies the test of heterogeneity did not 
reach statically significance (P- interaction = 0.06 in (Greer et al, 
2006) and P- interaction = 0.15 in our study). A limitation of the 
study by Greer et al (2006) was the inclusion of pre- and 
postmenopausal women and the lack of information on HT use 
among postmenopausal women. It has been reported that BMI and 
HT may interact on the risk of ovarian cancer (Leitzmann et al, 
2009; Collaborative Group on Epidemiological Studies of Ovarian 
Cancer, 2012). We were able to look at this interaction separately 
among nulliparous and parous women, and found some 
differences in the magnitude of the associations, although the 
number of cases in these subgroups was too small to reach firm 
conclusions. However, it may be worthwhile to assess this 
association in larger studies. 

In a case-control study the use of OC did not reduce the risk of 
ovarian cancer in nulliparous women, but it was associated with a 
lower risk in parous women (Ness et al, 2001). We did not observe 
any qualitative difference in the risk estimates for OC use between 
nulliparous and parous women. 

A previous prospective study based on 181 ovarian cancer cases 
found that the association between nulliparity and ovarian cancer 
was stronger among those with a first degree family history of 
breast or ovarian (Vachon et al, 2002). While first degree family 
history of breast or ovarian cancer was more strongly associated 
with ovarian cancer in nulliparous than in parous women in our 
study, neither association reached statistical significance and our 
estimates were lower than those found earlier (Vachon et al, 2002). 
The study by Vachon et al (2002) was also a prospective study of 
postmenopausal women with similar follow-up time, but our 
combined data set was larger (37 377 vs 125 437) and there were 
more incident epithelial ovarian cases (181 vs 623). However, there 
were a higher proportion of nulliparous women with missing 
information on family history than parous women in our study. 

Our study has several limitations. Despite combining two large 
studies, the number of cases was small and may have precluded us 
from observing statistically significant associations. However, our 
sample size was larger than any other previous study. Given the 
low incidence of ovarian cancer, a small number of cases is a 
common problem of prospective studies of this disease. Another 
limitation was that we did not have fertility information in the 
NIH-AARP study, and, therefore, we could not adjust for this 
variable, which may be particularly important in the nulliparous 
group. About 60% of the nulliparous women were ever married 
(married, widowed, divorced, or separated) or living as married. As 
these women might have reproductive and hormonal character- 
istics different from nulliparous women who were never married or 
did not live as married, we conducted sensitivity analyses restricted 
to women who were ever married or living as married. However, 
the results were similar to those presented in our tables. For some 
factors that have been associated with ovarian cancer, such as tubal 
ligation, endometriosis or infertility, information was only available 
in the PLCO data set, and the number of ovarian cancer cases 
among nulliparous women was too small (N= 8) to determine the 
ovarian cancer risk associated with them. Finally, we did not have 
follow-up information regarding oophorectomy after enrolment 
and, therefore, we could not censor women who removed their 
ovaries after the beginning of the study. Our study has several 
notable strengths such as the prospective design with long follow- 
up, a larger number of nulliparous women compared with previous 
studies, and the comprehensive assessment of hormonal risk 
factors. 

In summary, we did not observe significant differences in 
strengths of associations of known risk factors between among 
nulliparous and parous women. There was a suggestive interaction 
between BMI and parity in ever users of HT. Pooled analyses or 
larger studies are needed to provide more definitive and reliable 
evidence of this potential interaction. 
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Table 3. Risk of invasive epithelial ovarian cancer in relation to hormone-related factors by parity 



Nulliparous women Parous women 


Risk factor 


Cases 
(N-102) n(%) 


HR a 


(95% Cl) a 


Cases 
(N=521) n (%) 


HR b 


(95% Cl) b 


P-heterogeneity c 


Duration of use of OC 














0.98 


Never users or <1 year 
1-9 years 
^10 years 


76 (80.0) 
15 (15.8) 
4 (4.2) 


1.00 
0.78 
0.53 
P trend = 0.16 


ref. 
(0.44, 1 .39) 
(0.19, 1.47) 


333 (67.4) 
1 33 (26.9) 
28 (5.7) 


1.00 
0.81 
0.45 
P trend < 0.001 


ref. 
(0.66, 1 .00) 
(0.30, 0.66) 




Duration of HT use 














0.18 


Never users 
< 1 0 years 
^ 1 0 years 


55 (57.9) 

o a /or o\ 

24 (25.3) 
16 (16.8) 


1.00 
0.85 
1.46 

P,rend = 0.41 


ref. 
(0.52, 1 .40) 
(0.83, 2.58) 


1 95 (39.5) 

1 96 (39.7) 
103 (20.9) 


1.00 
1 .46 
1.61 

P trend < 0.001 


ref. 
(1.19, 1.79) 
(1.27, 2.06) 




Age at menarche 














0.24 


Early 

Middle 

Late 


46 (48.4) 
45 (47.4) 
4 (4.2) 


1.00 

1.10 
0.45 
P trend = 0.42 


ref. 
(0.73, 1 .6)] 
(0.16, 1.25) 


201 (40.7) 
232 (47.0) 
61 (12.4) 


1.00 
0.99 
0.85 
P trend = 0.38 


ref. 
(0.82, 1.21) 
(0.63, 1.14) 




Age at last menstrual cycle (years) 














0.48 


^44 
45-49 
50-54 
^55 


24 (25.3) 
29 (30.5) 
34 (35.8) 
8 (8.4) 


1.00 
1.13 
0.92 
1.24 

Ptrend = 0.96 


ref. 
(0.65, 1.94) 
(0.54, 1.56) 
(0.55, 2.78) 


137 (27.7) 
103 (20.9) 
202 (40.9) 
52 (10.5) 


1.00 
0.95 
1.15 
1.08 

Ptrend = 0.23 


ref. 
(0.73, 1.23) 
(0.92, 1 .44) 
(0.78, 1.50) 




Lifetime ovulatory cycles d (LOC) 














0.50 


Q4 (>470) 

Q3 (>41 1 and ^470) 
Q2 (>330 and ^411) 
Q1 (^330) 


37 (39.0) 
27 (28.4) 
25 (26.3) 
6 (6.3) 


1.00 

1.18 
1.10 
0.74 
P trend = 0.90 


ref. 
(0.72, 1.96) 
(0.65, 1.88) 
(0.29, 1.86) 


1 50 (30.4) 

1 1 7 (23.7) 

118 (23.9) 
109 (22.1) 


1.00 

0.81 
0.77 
0.76 
P trend = 0.05 


ref. 
(0.63, 1.03) 
(0.60, 1 .00) 
(0.58, 1.01) 




1 — j. 1 r - 1 1 - . r 1 . 

First degree family history of breast 
and/or ovarian cancer 














0.60 


No 
Yes 


48 (71.6) 
19 (28.4) 


1.00 
1.37 


ref. 
(0.81, 2.34) 


344 (76.6) 
105 (23.4) 


1.00 
1.21 


ref. 
(0.97, 1.51) 




BMI (kgm- 2 ) e 














0.15 


<25 

25-29.9 

^30 


39 (41.1) 
32 (33.7) 
24 (25.3) 


1.00 
1.36 
1 .48 

P,rend = 0.32 


ref. 
(0.79, 2.36) 
(U.o I , o.bo] 


237 (48.0) 
144 (29.2) 

110 m o\ 
I I 3 (22. V) 


1.00 

0.84 
1 .03 

Ptrend = 0.71 


ref. 
(0.66, 1.07) 
(U.6V, \ .06} 




Hysterectomy 














0.54 


No 
Yes 


82 (86.3) 
13 (13.7) 


1.00 

0.83 


ref. 
(0.45, 1.53) 


346 (70.0) 
148 (30.0) 


1.00 
1.00 


ref. 
(0.82, 1.23) 




Smoking status 














0.10 


Never 

Former 

Current 


43 (45.7) 
35 (37.2) 
16 (17.0) 


1.00 
0.99 
1.40 


ref. 
(0.63, 1.55) 
(0.78, 2.49) 


253 (52.0) 
189 (38.8) 
45 (9.2) 


1.00 
0.93 
0.76 


ref. 
(0.77, 1.13) 
(0.55, 1.05) 




Abbreviations: BMI = body mass index; CI = confidence interval; HR = hazard ratio; HT= hormone therapy; OC = oral contraceptive. Numbers only for women with no-missing values for the 
exposure of interest and corresponding adjusting variables. Age at menarche groups for the AARP data set: Early, ^ 12 years; Middle, 13-14 years; Late, ^15 years. Age at menarche groups for 
the PLCO data set: Early, <11 years; Middle 12-13 years; Late, >14 years. 

a Adjusted for BMI (continuous), duration of use of oral contraceptives (never or < 1/1-9/^ 10 years) and duration of use of hormone therapy (Never/< 10/^ 10 years), first degree family history 
of breast and/or ovarian cancer (no/yes/missing) and stratifying the baseline function by study (PLCO vs AARP). 
^Adjusted for same variables as nulliparous analysis plus number of live births (1/2/3-4/^5). 

c Based on a likelihood-ratio test comparing a model fitted to all women with variables as in the parous analysis plus parity status (yes/no) and an interaction term of parity status time the 
exposure of interested (coded in continuous form) to a model without the interaction term. 

^Ovulatory cycles computed according to the formula [(age at menopause — age at menarche-time spent pregnant-duration of OC use)*365.25]/average cycle length, as given in Cramer et al 
(1995). The average cycle length was set to 28 days. All the other characteristics were set to the mid-point of the interval of the corresponding variable, separately for the PLCO and for the NIH- 
AARP. Quartiles of lifetime ovulatory cycles were defined based on the distribution in the entire study population. 
e At enrolment. 
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Table 4. Risk of invasive epithelial ovarian cancer in relation to BMI by HT use and parity 



Nulliparous women 



Parous women 



Stratum/risk factor Cases n (%) 
Never users of HT 

BMI (kgm- 2 ) 



(95% Cl) a 



Cases n (%) 



HR 



(95% CI)" 



P-heterogeneity c 



0.76 



<25 

25-29.9 

>30 



Ever users of HT 



20 (36.4) 

21 (38.2) 
14 (25.5) 



1.00 
1.90 
2.10 



ref. 
(0.86, 4.19) 
(0.53, 8.27) 
P, rend = 0.20 



70 (35.9) 
59 (30.3) 
66 (33.9) 



1.00 
1.10 
1.83 



ref. 
(0.73, 1.66) 
(0.95, 3.52) 

P,re„d = 0.11 



BMI (kgm"1 

<25 
25-29.9 

5=30 



19 (47.5) 
11 (27.5) 
10 (25.0) 



1.00 
1.01 
1.50 



ref. 
(0.47, 2.18) 
(0.57, 3.92) 
P trend = 0.548 



167 (55.9) 
85 (28.4) 
47 (15.7) 



1.00 
0.74 
0.70 



ref. 
(0.55, 0.99) 
(0.44, 1.11) 
P trend = 0.05 



0.08 



Abbreviations: BMI = body mass index; C! = confidence interval; HR = hazard ratio; HT= hormone therapy; PLCO = Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial. Numbers 
only for women with no-missing values for the exposure of interest and corresponding adjusting variables. 

a Adjusted for BMI (continuous), duration of use of oral contraceptives (never or < 1/1-9/^ 10 years), first degree family history of breast and/or ovarian cancer (no/yes/missing) and stratifying 
the baseline function by study (PLCO vs AARP). 

^Adjusted for same variables as nulliparous analysis plus number of live births (1/2/3-4/^5), age at first birth (< 16/1 6-1 9/20-24/25-29/30-34/^ 35 years} and duration of use of hormone 
therapy (Never/<10/3s10 years). 

c Based on a likelihood-ratio test comparing a model fitted to all women with variables as in the parous analysis plus parity status (yes/no) and an interaction term of parity status times the 
exposure of interested (coded in continuous form} to a model without the interaction term. 
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